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O \ Abstract 

We discuss the possibilities of high precision measurement of the solar neutrino mixing 
JL , angle 9q = 9\2 in solar and reactor neutrino experiments. The improvements in the de- 

9 

termination of sin 6\2, which can be achieved with the expected increase of statistics and 
reduction of systematic errors in the currently operating solar and KamLAND experiments, 
are summarised. The potential of LowNu v — e elastic scattering experiment, designed to 
measure the pp solar neutrino flux, for high precision determination of sin 2 9±2, is investigated 
in detail. The accuracy in the measurement of sin 2 #12, which can be achieved in a reactor 
experiment with a baseline L ~ (50 — 70) km, corresponding to a Survival Probability MIN- 
imum (SPMIN), is thoroughly studied. We include the effect of the uncertainty in the value 
of sin 2 #13 in the analyses. A LowNu measurement of the pp neutrino flux with a 1% error 
would allow to determine sin 2 612 with an error of 14% (17%) at 3a from a two-generation 
(three-generation) analysis. The same parameter sin 2 612 can be measured with an uncer- 
tainty of 2% (6%) at la (3a) in a reactor experiment with L ~ 60 km, statistics of ~60 
GWkTy and systematic error of 2%. For the same statistics, the increase of the systematic 
error from 2% to 5% leads to an increase in the uncertainty in sin 2 812 from 6% to 9% at 
3a. The inclusion of the sin 2 #13 uncertainty in the analysis changes the error on sin 2 612 
to 3% (9%). The effect of sin 2 #13 uncertainty on the sin 2 612 measurement in both types of 
experiments is considerably smaller than naively expected. 
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1 Introduction 



There has been a remarkable progress in the studies of neutrino oscillations in the last several years. 
The experiments with solar, atmospheric and reactor neutrinos [1, 2, 3, 4, 5, 6, 7] have provided 
compelling evidences for the existence of neutrino oscillations driven by non-zero neutrino masses 
and neutrino mixing. Evidences for oscillations of neutrinos were obtained also in the first long 
baseline accelerator neutrino experiment K2K [8] . 

The recent Super-Kamiokande data on the L/i?-dependence of multi-GeV //-like atmospheric 
neutrino events [6], L and E being the distance traveled by neutrinos and the neutrino energy, and 
the new more precise spectrum data of KamLAND and K2K experiments [9, 10], are the latest 
significant contributions to this progress. For the first time the data exhibit directly the effects 
of the oscillatory dependence on L/E and E of the probabilities of //-oscillations in vacuum [11]. 
We begin to "see" the oscillations of neutrinos. As a result of these magnificent developments, 
the oscillations of solar is e , atmospheric and 9^, accelerator (at L ~250 km) and reactor 
v e (at L ~180 km), driven by nonzero //-masses and //-mixing, can be considered as practically 
established. 

The SK atmospheric neutrino and K2K data are best described in terms of dominant 2-neutrino 
Vn — > v T (Pfj, — > u T ) vacuum oscillations. The best fit values and the 99.73% C.L. allowed ranges of 
the atmospheric neutrino oscillation parameters |Am^| = jAm^l and 9 a = #23 read [6]: jAmiy = 
2.1 x 10" 3 eV 2 , sin 2 2^ 2 3 = 1.0, [Am^l = (1.3 - 4.2) x 10" 3 eV 2 , sin 2 29 23 > 0.85. The sign of 
Am 2 x and of cos 26> 2 3, if sin 2 29 2 3 ^ 1-0, cannot be determined using the existing data. 

The combined 2-neutrino oscillation analysis of the solar neutrino and the new KamLAND 766.3 
Ty spectrum data shows [9, 12, 13] that the solar neutrino oscillation parameters lie in the low- 
LMA region : Am| = Am 2 21 = (7.9±° 6 5 ) x 10" 5 eV 2 , tan 2 9 Q = tan 2 9 12 = (0.401&0?) [9]. The 
high-LMA solution is excluded at more than 3a. The value of Am^ is determined with a remark- 
ably high precision of 12% at 3a. Maximal solar neutrino mixing is ruled out at ~ 6a. 

The solar and atmospheric neutrino, and KamLAND and K2K neutrino oscillation data require 
the existence of three-neutrino mixing in the weak charged lepton current. In this case the neutrino 
mixing is characterised by one additional mixing angle #13 - the only small mixing angle in the 
PMNS matrix. Three-neutrino oscillation analyses of the solar, atmospheric and reactor neutrino 
data show that sin 2 13 < 0.05 [12, 13] l . 

Understanding the origin of the patterns of solar and atmospheric neutrino mixing and of 
Ara^ and Am| l5 suggested by the data, is one of the central problems in neutrino physics today. 
A pre-requisite for any progress in our understanding of neutrino mixing is the knowledge of the 
precise values of solar and atmospheric neutrino oscillation parameters, #12, Am^, and $23, Am^ 
and of #i3 . In the present article we discuss the possibilities of high precision measurement of the 
solar neutrino mixing angle 612 in solar and reactor neutrino experiments. 

The solar neutrino mixing parameter sin 2 9 12 is determined by the current KamLAND and solar 
neutrino data with a relatively large uncertainty of 24% at 3a. In the future, more precise spectrum 
data from the KamLAND experiment can lead to even more accurate determination of the value 

1 After the new background data published by the KamLAND collaboration [9] the bound changes to sin 2 13 < 
0.055 [12] whereas inclusion of the recently published 391 days SNO salt phase data [14] gives sin 2 #13 < 0.044 [15] 
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of Am 2 ^. However, these data will not provide a considerably more precise measurement of sin 2 9 12 
owing to the fact that the baseline of the KamLAND experiment effectively corresponds to a v e 
Survival Probability MAXimum (SPMAX) [16, 17]. The analysis of the global solar neutrino data 
taking into account a possible reduction of the errors in the data from the phase-Ill of the SNO 
experiment shows that the uncertainty in the value of sin 2 9 12 would still remain well above 15% 
at 3a [18]. 

We begin by summarising the results on sin 2 6^2, obtained using the current global solar and 
reactor neutrino oscillation data (section 2). We consider the improvements in the determina- 
tion of sin 2 # 12 , which can be achieved with the expected increase of statistics and reduction of 
systematic errors in the experiments which are currently operating. In particular, the effect of 
KamLAND data, corresponding to a statistics of 3 kTy, as well as of the data of phase-Ill of SNO 
experiment, are analysed. 

We turn next to future experiments. We discuss first (in section 3) the possibility of high 
precision determination of sin 2 9u in a LowNu solar neutrino experiment, designed to measure the 
~MeV and sub- Me V components of the solar neutrino flux: pp, pep, CNO, 7 Be. It is usually 
suggested that the LowNu experiments can provide one of the most precise measurements of the 
solar neutrino mixing angle [19, 20, 21]. Detailed analysis was carried out in [22] and it was 
concluded that a future pp experiment should have accuracy better than 3% in order to improve 
on the knowledge of tan 2 12 . 

We consider a generic v — e scattering experiment measuring the pp neutrino flux and perform 
a detailed quantitative analysis of the precision with which sin 2 9yi can be determined in such an 
experiment. We examine the effect of including different representative values of the pp neutrino 
induced event rate in the x 2 analysis of the global solar neutrino data. Three values (0.68, 0.72, 
0.77) of the (normalized) event rate from the currently allowed 3a range are considered. The error 
in the measured rate is varied from 1% to 5%. We investigate how much the accuracy on sin 2 9\2 
improves by adding the pp flux data in the analysis. The dependence of the sensitivity to sin 2 9 12 
on the central value of the measured pp flux as well as on the measurement errors is studied. 
We compare the precision in sin 2 9yi expected with assumed data on pp neutrinos included in 
the analysis, with the sensitivity that can be achieved using prospective results from phase-Ill of 
SNO experiment and 3 kTy KamLAND data, and comment on the minimum error required in 
the LowNu pp experiments to improve the precision of sin 2 9yi measurement. The impact of the 
uncertainty due to #13 on the allowed ranges of sin 2 9yi is studied as well. 

In the section 4 we analyse in detail the possibility of a high precision determination of sin 2 9\2 
in a reactor experiment with a baseline of L ~ (50 — 70) km, corresponding to a v e Survival 
Probability MINimum (SPMIN). That such an experiment can provide the highest precision in 
the measurement of sin 2 6*12 was pointed out first in [16]. A rather detailed study of the preci- 
sion in sin 2 # 12 , which might be achieved in an SPMIN experiment with the flux of v e from the 
Kashiwazaki reactor complex in Japan and L = 54 km, was performed recently in [23] 2 . We 
consider a generic SPMIN reactor experiment with a KamLAND-type detector. We investigate 
the dependence of the precision on sin 2 #12 which can be achieved in such an experiment on the 
baseline, statistics and systematic errors. More specifically, the spectrum data is simulated for 

2 In ref. [23] this experiment is called SADO. 
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four different true values of Am^ and for each of these values the optimal baseline at which the 
most precise measurement of sin 2 12 could be performed is determined. We show, in particular, 
that an independent determination of Ara^ with sufficiently high accuracy would allow sin 2 9 i2 
to be measured with the highest precision over a relatively wide range of baselines. The effect of 
sin 2 9 13 uncertainty on the sin 2 9 12 determination is investigated in detail. 
The results of the present study are summarised in section 5. 

2 Measuring sin 2 612 in Existing Experiments 

In this section we review the precision of Am 2 x and sin 2 9\ 2 determination from the existing solar 
neutrino and KamLAND data from two-neutrino oscillation analysis. We also discuss possible 
improvements in the precision that can be achieved in the currently running experiments. 

2.1 Current Solar and Reactor Neutrino Data on sin 2 #12 
In the present global solar neutrino and KamLAND data analysis we include 

• data from the radiochemical experiments, CI [1] and Ga (Gallex, SAGE and GNO combined) 
[2], 

• the 1496 day 44 bin Zenith angle spectrum data from SK [3], 

• the 34 bin combined CC, NC and Electron Scattering (ES) energy spectrum data from the 
phase I ( pure D 2 phase) of SNO [4], 

• the data on CC, NC and ES total observed rates from the phase II (salt phase) of SNO 
experiment [5]. 

The 8 B flux normalization factor is left to vary freely in the analysis, while for the pp, pep, 
7 Be, CNO, and hep fluxes the predictions and uncertainties from the recent standard solar model 
(SSM) [24] (BP04) are used. We skip the details of the x 2 ~ analysis, which can be found in 
[25, 26, 27]. 

We include the 766.3 Ty KamLAND data in the global analysis. For treatment of the latest 
KamLAND data in the combined analysis we refer the reader to [12], while details regarding the 
future-projected analysis with simulated KamLAND data are given in [28]. The best-fit in the 
combined analysis of solar and KamLAND data is obtained for [12] 3 

• Am 2 x = 8.4 x 1(T 5 eV 2 , sin 2 9 12 = 0.28, f B = 0.88 

In Fig. 1 we plot A% 2 = X 2 ~Xmin as a function of Am^ (right panel) and sin 2 9 12 (left panel), 
for a two-neutrino oscillation fit of the global solar neutrino + KamLAND data. The parameters 

3 With the inclusion of the new KamLAND background data [9], the best- fit shifts to slightly smaller values of 
Am^: Am^ = 8.0 x 10" 5 cV 2 , sin 2 6 12 = 0.28, f B = 0.88 [12]. The inclusion of the recent SNO results give 
the best-fit parameters as Am^ = 8.04 x 10~ 5 cV 2 , sin 2 12 = 0.31, [15] 
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Figure 1: A% 2 as a function of sin 2 9 12 (left panel) and of Am| x (right panel). The curves shown are 
obtained from global analysis of the current solar neutrino data (black dot-dashed line), current 
solar + KamLAND data (red solid line), solar neutrino data with projected SNO-III errors + 
current KamLAND data (green short dashed line) and the solar neutrino data with projected 
SNO-III errors + prospective 3 kTy KamLAND data (blue long dashed line). The horizontal lines 
indicate the 3a limit (A% 2 = 9) for 1 parameter fit. 



not given in each of the two panels are allowed to vary freely. From this figure one can easily read 
off the allowed ranges of the displayed parameter at various confidence levels. The horizontal line 
shows the 3a limit corresponding to A% 2 = 9 for a 1 parameter fit. In Table 1 we present the 
current 3a allowed ranges of Am^ and sin 2 9 12 obtained from Fig. 1. We also give the spread in 
these parameters which is defined as 

spread = x 100 (1) 

where prm is either Am^i or sin 2 9u- Table 1 and Fig. 1 demonstrate clearly that the KamLAND 
experiment has remarkable sensitivity to Amf j . The inclusion of KamLAND data with increasing 
statistics in the analysis progressively reduces the spread in Am^ , and with the latest data the 
3a spread is ~ 12%. This demonstrates, in particular, the extraordinary precision that has already 
been achieved in the determination of Am^. On the other hand, the spread in sin 2 9 12 is seen 
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to be controlled mainly by the solar neutrino data and does not show any marked reduction 
with the inclusion of the KamLAND data in the analysis. Using the global solar neutrino and 
KamLAND data allows to determine sin 2 9 U with an error of 24% at 3a. 



2.2 Prospective KamLAND and SNO Data and sin 2 #i 2 

We will analyze next the expected impact of future data from KamLAND and SNO experiments 
on the Aula and sin 2 Qyi determination. 

The SNO experiment is sensitive to the flux of 8 B neutrinos with energy E ^ 5 MeV. For the 
oscillation parameters in the low-LMA region, the 8 B neutrino survival probability of interest is 
given approximately by (MSW adiabatic transition probability): 



The CC/NC ratio measured in SNO determines P ee independently of the 8 B neutrino flux nor- 
malization. Thus, it can give a direct measure of sin 2 6*12. Reducing the errors (A(R C c/Rnc) 
in the measured CC and NC event rates in SNO, Rcc an d Rnc, can improve the precision of 
determination of this parameter since 



The oscillations of reactor v e detected in KamLAND experiment are practically not affected 
by Earth matter effects and the corresponding v e survival probability has the form 



The average energy and baseline for KamLAND correspond to shi^Am 2 ^ L/AE) ps 0, i.e., to 
a Survival Probability MAXimum (SPMAX). As a consequence, the coefficient of the sin 2 26 12 
term in P^ L is relatively small, weakening the sensitivity of KamLAND to Ou- As was shown in 
[16], the most precise measurement of sin 2 # 12 can be performed in a reactor experiment with a 
baseline tuned to a Survival Probability MINimum (SPMIN), i.e., to shi^Am^ L/AE) as 1. We 
will discuss the sensitivity to sin 2 9 U , which can be achieved in such an experiment, in section 4. 

In phase-Ill of the SNO experiment, the NC events will be observed directly (and independently 
from the CC events) using 3 He proportional counters. This will help to increase the NC statistics 
and reduce the systematic errors in the NC data. In addition, the correlations between the errors 
in the measured CC and NC event rates will be absent. The total projected error in the measured 
NC event rate in phase-Ill of SNO experiment is ~ 6% [29]. We incorporate this in our analysis 
instead of the present error in R NC of 9%. For the CC event rate Rcc measured at SNO we 
assume a somewhat reduced total error of 5% (the current error in Rcc is approximately 6% 
[14]). We assume also that the central values of the measured CC and NC even rates will remain 
unchanged. The results of our analysis are shown in Fig. 1. The short dashed lines in Fig. 1 
display the behavior of A% 2 with anticipated SNO phase-Ill results and the reduced projected 
errors added to the global solar neutrino and KamLAND 766.3 Ty spectrum data. The figure 



P ee ( 8 B )^sm 2 9 12 . 



(2) 



A(sin 2 fl 12 ) = A(R CC /R NC ). 



(3) 




(4) 
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shows that with the inclusion of the SNO phase-Ill prospective results, the allowed range of Am^ 
in the low-LMA region remains unchanged: it is determined principally by the KamLAND data. 
However, the higher Am^ regions get more disfavored as the reduced errors in Rcc and Rnc lead 
to a stronger rejection of larger values of Rcc /Rnc [30]. The figure also shows that the allowed 
range of sin 2 12 gets further constrained from above. From Table 1 we see that the 3a spread of 
sin 2 6*i2 becomes 21% with the inclusion of projected SNO phase-Ill data in the presently existing 
set of data. 

We also studied the effect of increased statistics of KamLAND experiment on the sin 2 # 12 and 
Am 2 ! determination. To this end, we include 3 kTy KamLAND spectrum data, simulated at 
Am 21 = 8.3 x 10~ 5 eV 2 , in our analysis with projected SNO phase-Ill data. We use a systematic 
error of 5% for KamLAND since the KamLAND systematic error is expected to diminish 4 after 
the planned fiducial volume calibration and re-evaluation of the uncertainties in the power of 
the relevant nuclear stations. If the real KamLAND spectrum data conforms to this projected 
spectrum, the allowed range of Am 21 would be further constrained, allowing the determination of 
Am 21 with an accuracy of about 5%. The higher Am 21 regions would be disfavored even stronger. 
The uncertainty in the value of sin 2 9 U would be smaller and the 3a spread, as seen from Table 
1, could be 18%. Clearly, the uncertainty in sin 2 9 12 cannot be reduced to 15% or less (at 3a) by 
future data from the currently operating solar and reactor neutrino experiments 5 



Data set 


(3a) Range of 


(3a) spread in 


(3a) Range of 


(3a) spread in 


used 


Ara^j eV 2 




sin 2 


#12 


sin 2 Q 12 


only sol 


3.3 - 15.3 


65% 


0.22- 


■0.38 


27% 


sol+ 766.3 Ty KL 


7.4 - 9.5 


12% 


0.22- 


■0.36 


24% 


sol(+SN03) + 766.3 Ty KL 


7.4 -9.5 


12% 


0.22 - 


0.34 


21% 


sol(+SN03)+3KTy KL 


7.7 -8.9 


7% 


0.23 - 


0.33 


18% 



Table 1: The 3a allowed ranges and spread (in per cent) of Am 21 and sin 2 6 l 12 obtained from 1 
parameter fits. 



3 Determining 612 from Measurement of the pp Solar Neu- 
trino Flux 

The pp fusion reaction is the main contributor to the observed solar luminosity and the correspond- 
ing pp neutrinos constitute the largest and dominant component of the solar neutrino flux. The 
SSM uncertainty in the predicted solar neutrino fluxes is the least for the pp neutrinos (~1%) [24]. 
So far only the Ga experiments have provided information on the pp neutrino flux because of their 

4 The KamLAND systematic error at present is 7.13% including the (a-n) background and 6.5% without this 
background [9]. 

5 If we include the recently published SNO results [14] in our analysis then the 3a spread in sin 2 6>i 2 be- 
comes 26% whereas the solar+KamLAND analysis including the recent SNO results and the (a, n) background 
in KamLAND [9] gives the spread as 22% [15]. 
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low energy threshold of 0.23 MeV. Sub-MeV solar neutrino experiments (LowNu experiments) 
are being planned for measuring the flux of pp neutrinos using either charged current reactions 
(LENS, MOON, SIREN) or v - e~ elastic scattering (XMASS, CLEAN, HERON, MUNU, GE- 
NIUS) [19, 20, 21]. Since according to the SSM, most of the energy released by the Sun (~99%) 
is generated in the pp— cycle of reactions in which also the pp neutrinos are produced, a precise 
determination of the pp neutrino flux would lead to a better understanding of the solar energetics 
and, more generally, of the physics of the Sun. It has also been realized that high precision mea- 
surement of the pp neutrino flux can be instrumental for more accurate determination of the solar 
neutrino mixing parameter, which, as we have seen in the preceding section, will not be determined 
with an uncertainty smaller than ~18% (at 3a) by the currently operating experiments. 

Since the pp neutrino energy spectrum extends upto 0.42 MeV only, for Am^ in the LMA 
region, the pp neutrino oscillations are practically not affected by matter effects in the Sun or 
the Earth. Thus, to a good approximation, the pp neutrino oscillations are described by the v e 
survival probability in the case of oscillations in vacuum, in which the oscillating term is strongly 
suppressed by the averaging over the region of neutrino production in the Sun [31]: 

Pf;(H-l-isin 2 2£ 12 . (5) 

The normalised event rate for a v — e scattering (ES) and charged current (CC) experiments 
measuring the pp neutrino flux is given respectively by 

R pp = i P ee(PP)) + r ppi l ~ ( P ee(PP))) for ES experiments, (6) 
R pp = (PeeiPP)) f° r CC experiments, , (7) 

where r pp = J f pp (E)a{v IJt i T \e~)dE/ J f pp (E)a(v e e~)dE m 0.3 where f pp {E) is the pp-neutrino 
flux, cr(z/ M ( r )e~) and <r{v e e~) are the v^t) — e ~ and u e — e~ eleastic scattering cross sections and 
(..) denotes averaged probabilites. 6 

For an ES experiment, the second term in eq. (6) represents the NC contribution. Since the pp 
neutrino survival probability is largely independent of energy one can use eq. 5 for the averaged 
probabilites . Using eqs. (5) - (7), one finds for the uncertainty in sin 2 dyi determination: 

A R 

A(sin 2 M^f {£- for CC experiments (8) 

pp 2 cos 26*12 

AD 1 

A(sin 2 12 )gf ~ - ^ for ES experiments (9) 

2 COS 2c/i2 1 — T pp 

where AR pp is the error in the measured value of R pp . A comparison of eq. (8) with eq. (3) 
shows that for the same value of AR pp and of the error in the CC to NC event rate ratio measured 
by SNO, A(Rcc I Rnc) , the pp neutrino experiments of the CC type can provide a more precise 
measurement of sin 2 9 12 only if cos 29^ > 0.5 (sin 2 #12 < 0.25). Similarly, it follows from eq. (9) 
and eq. (3) that for AR PP = A(Rcc / Rnc) , the ES pp experiment could provide a more precise 

6 Strictly speaking, eq. (6) is an approximate expression since P^(pp) depends on the neutrino energy, while 
eq. (6) is valid for energy-independent P^(pp). In our numerical analysis we have not used this approximation. . 
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determination of sin 2 9n only if the true value of cos2#i 2 > 0.71 (sin 2 9\ 2 < 0.14). Since the 
currently allowed 3a range of sin 2 9 12 is sin 2 12 = (0.22 — 0.38), for almost all of the allowed values 
of sin 2 #i 2 , SNO will have a better sensitivity to sin 2 12 than a LowNu pp experiment measuring 
the pp neutrino flux with the same experimental error as the error in the SNO data on the ratio 
Reel Rnc- 

In order to improve the accuracy of sin 2 9 12 determination after the SNO phase-Ill results will 
be available, the total experimental error in the measured event rate in the pp experiments has 
to be sufficiently small, which requires high statistics and well understood systematics. In the 
present section we will quantify these statements by incorporating hypothetical data on the pp 
neutrino flux in our analysis. We consider a generic pp neutrino ES experiment 7 and consider 
some illustrative sample rates from the currently predicted range. We give quantitative estimate 
of the sensitivity to sin 2 #i 2 expected to be achieved in i) a pp experiment, and ii) combining 
the prospective data from a pp experiment with the global solar and reactor neutrino data. In 
particular, we estimate the maximal error in a pp flux measurement, for which the uncertainty in 
the determined value of sin 2 9\2 would be smaller than that expected after the inclusion of the SNO 
phase-Ill results. We use in this analysis the pp neutrino flux and its 1% uncertainty predicted by 
the BP04 SSM [24]. The 1% error due to SSM uncertainties is added to the experimental errors 
in the pp flux determination. 

3.1 Two Generation Analysis 

We consider a generic v — e scattering experiment which can measure the pp neutrino flux. The 
experiment is assumed to have e~ kinetic energy threshold of 50 keV. We suppose that the BP04 
SSM predicts correctly the pp flux. The predicted event rate ( u pp rate") in such an experiment 
for the best fit values of the oscillation parameters, normalized to the rate predicted by the BP04 
model in the absence of neutrino oscillations, is 0.71; the predicted 3a range for the normalized 
pp rate is 0.67 — 0.76 [12]. We will consider three illustrative values of the normalized pp rate, 
0.68, 0.72 and 0.77, and vary the experimental error in the measured pp rate from 1% to 5%. The 
estimated theoretical uncertainties due to the SSM and their correlations for the pp neutrino flux 
are included in the analysis following the standard covariance approach [32]. Thus, we minimise 
the x 2 defined as 

N 

xi = E (^ cxpt - iC cory )(4rW pt - R t cmy ) ( 10 ) 

where Ri are the solar neutrino data points, N is the number of data points and (of,) -1 is the 
inverse of the covariance matrix, containing the squares of the correlated and uncorrelated exper- 
imental and theoretical errors. The 8 B flux normalisation factor f B is left to vary freely in the 

comparison of eqs. (8) and (9) shows that the CC pp experiments could achieve a better sensitivity on 
sin 2 6*12 due to the absence of the « 30% "contamination" caused by the NC events present in the ES event sample. 
In what follows, we present results for an ES pp neutrino experiment (XMASS, etc.). Henceforth the term "pp" 
experiment implies an ES pp neutrino experiment. 
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Figure 2: A% 2 as a function of sin 2 9 12 . The results shown are obtained in a two-neutrino oscillation 
analysis of the global solar neutrino and KamLAND 766.3 Ty spectrum data and simulated data 
from a LowNu pp neutrino experiment. The three panels correspond to three illustrative values of 
the event rate due to pp neutrinos in the LowNu experiment ("pp rate"), normalized to the rate 
predicted by the BP04 SSM. In each case results for four different assumed values of the error 
in the measured pp rate are shown. We also show the curve obtained in global analysis of the 
current solar and reactor neutrino data. The horizontal line indicates the 3a limit (A% 2 = 9) for 
1 parameter fit. 

analysis. The errors and correlations due to the other fluxes are taken from the BP04 SSM 8 . 

Let us begin by analysing first the potential results from a possible future pp neutrino experi- 
ment alone. In this case eq. (6) can be used to get the approximate values of sin 2 29i 2 for a given 
pp rate, 

sm 2 26 12 ^ 2{1 ~ Rpp \ (11) 

If we assume a pp rate of 0.72 and 1% experimental error (in addition to the SSM uncertainty of 
1%), the x 2 analysis gives for the 3a range of allowed values of sin 2 9\ 2 0.21 < sin 2 9\ 2 < 0.33. This 
agrees very well with what one would obtain using eq. (11) including the errors along with the 
mean pp rate value of 0.72. Thus, the spread in sin 2 9± 2 is of about 22%, which is not much smaller 

8 One can also keep the normalization of the pp neutrino flux as a variable parameter, subject to the solar 
luminosity constraint [22]. 
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Figure 3: The la (68.27% C.L.), 1.64a (90% C.L.), 2a (95.45% C.L.) and 3a (99.73% C.L.) allowed 
range of sin 0\2 21 

continuous function of the error in the pp rate for three assumed mean values 
of the measured pp rate. Shown is also the current 3a allowed range of sin 2 9\2- 

than the spread in sin 2 # 12 determined using the current data. In what follows, we will perform an 
analysis of the data from all experiments, including illustrative rates from a v — e scattering type 
pp experiment , when estimating the sensitivity to the mixing angle 9yi- 

Figures 2 and 3 show results obtained in a two-neutrino oscillation analysis of the KamLAND and 
global solar neutrino data, including the pp rate assumed to be measured in the LowNu ES ex- 
periment. In Fig. 2 we plot the dependence of A% 2 on sin 2 # 12 and show results for four values 
of the error in the measured pp rate, 1%, 2%, 3% and 4%. The cyan colored lines correspond 
to results obtained using the currently existing data. One can easily read from the figure the 
range of allowed values of sin 2 #12 at a given C.L. for a chosen experimental error in the pp rate 
measurement. In Fig. 3 the corresponding allowed range of sin 2 # 12 is shown as a function of the 
error in the pp rate measurement. We let the pp rate error vary from 1% to 5%. The various bands 
correspond to 68.27% C.L. (la), 90% C.L. (1.64a), 95.45% C.L. (2a) and 99.73% C.L. (3a). For 
comparison, the 3a range of the presently allowed values of sin 2 6*12 is indicated in the figure by 
horizontal lines. 

In both figures we present results for the 3 assumed values of the measured pp rate, R pp . The 
figures show that: 

1. For R pp = 0.68 at the lower end of the predicted range of values of the pp rate, the lower 
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bound on sin 2 9u increases considerably with reducing the error in the pp flux measurement. 
The upper bound also increases but not significantly: the reduction of the error in pp rate 
measurement has a much smaller effect on the maximal allowed value of sin 2 9 12 . 

2. For R pp = 0.72, which is close to the best-fit predicted rate, both the lower bound on sin 2 On 
increases and the upper bound decreases, tightening the allowed range of sin 2 9\ 2 . Reducing 
the error in the pp flux measurement improves the precision of determination of sin 2 6*12. 

3. For the relatively high pp rate, R pp = 0.77, the minimal allowed value of sin 2 9 12 diminishes 
somewhat, while the maximal allowed value diminishes considerably. 

These features can be understood by analyzing the expression for the probability of survival of 
pp neutrinos given in eq. (6). It is evident that a lower pp rate drives sin 2 12 towards higher 
values and vice versa. Therefore the lower bound on sin 2 9 12 increases for R pp = 0.68. The effect 
of reducing the error in the pp flux measurement has the effect of pushing 9 12 towards higher 
values. Likewise, one could expect that the maximal value of sin 2 9i 2 should equally increase for 
R pp = 0.68. However, since we have used the global solar neutrino data in the analysis, the 
maximal allowed value of sin 2 9± 2 increases only slightly as the corresponding higher values of 
sin 2 9 12 are strongly disfavored by the already existing data. Thus, if a future pp (ES) experiment 
measures a value of R pp near the lower end of the presently predicted range, the lower limit on 
sin 2 9 12 will increase as the error in the measured R pp is reduced. The upper limit will increase 
slightly, but the effect of reducing the pp rate error will not be drastic. 

A higher measured value of R pp , R pp = 0.77, requires a lower value of sin 2 # 12 . Consequently, 
the maximal allowed value of sin 2 6 12 is seen to diminish substantially. The lower limit on sin 2 12 
could be pushed to relatively small values by the data from the pp experiment alone. However, 
such small values of sin 2 9\ 2 are already excluded by the current set of data and therefore the lower 
limit on sin 2 9\ 2 cannot reduce much. 

A pp rate of R pp = 0.72 is quite consistent with the current best-fit values of the parameters. 
As a consequence, the corresponding maximal allowed value of sin 2 9\ 2 diminishes and the mini- 
mal value increases as the error in R pp is reduced. Thus, the precision of sin 2 9\ 2 determination 
increases. 

We summarize in Table 2 (columns 3 and 4) the range of allowed values of sin 2 9 12 expected from 
a combined analysis of KamLAND and global solar neutrino data and the future (hypothetical) 
data from the pp experiment. With a 1% experimental error in the pp rate, the 3<r spread can 
decrease to about 14%. We note that the maximal and/or minimal allowed values of sin 2 9 12 
depend critically on the measured mean value of the pp rate. However, the spread in sin 2 9\2 
is practically independent of the latter. For R pp = 0.72, the 3<r spread in the value of sin 2 9 12 
is approximately 23%, 19%, 18% and 14% respectively for an error in the pp rate of 4%, 3%, 
2% and 1% (Table 2, 4th column). The 3a spread in sin 2 9 12 without including the phase-Ill 
SNO data is 24% (Table 1). We observe that pp flux measurement can improve the precision of 
sin 2 9 12 determination even with an error > 4-5%. In Fig. 4 we show the expected range of 
allowed values of sin 2 9 12 as a function of the error in the measured pp rate, after the potential 
results and projected errors from phase-Ill of the SNO experiment (SN03) are included in the 
global analysis. If the measured mean pp rate is 0.68 (0.77), the maximal (minimal) allowed value 
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Figure 4: Same as Fig. 3 but including the prospective data (and errors) from the phase-Ill of 
SNO experiment. 



of sin 2 9i2 slightly increases (decreases) and the minimal (maximal) value increases (decreases) 
significantly. The total uncertainty in sin 2 # 12 is independent of the central value of the pp rate, in 
agreement, with what we have found earlier. For R pp = 0.72, the 3cr spread in sin 2 9 12 is ~ 21%, 
18%, 18%, 14% for a 4%, 3%, 2% and 1% experimental error in R pp , respectively. 

It follows from our analysis that the data from the pp experiment can allow to reduce the 
expected uncertainty of 21% in the determination of sin 2 6 I 12 after the prospective phase-Ill SNO 
results are included in the analysis only if the error in the measured pp rate does not exceed ~ 
4%. However, even with ARp P = 1%, the 3a spread in the value of sin 2 9vz would not be smaller 
than ~14%. 

Table 2 summarises the results on the allowed ranges and spread of sin 2 #12, obtained in a 
two-neutrino oscillation analysis including SNO-III projected errors. As expected, the spread in 
sin 2 9 12 reduces with inclusion of the phase-Ill SNO data in the analysis. Note that the ranges 
given in Table 2 are with the present KamLAND 766.3 Ty spectrum data. If we use future higher 
statistics data from KamLAND , the allowed spread in sin 2 9 12 may reduce somewhat. 

3.2 The Impact of Non-Zero 13 

The solar and atmospheric neutrino, K2K and KamLAND data suggest the existence of 3-neutrino 
mixing and oscillations (see, e.g., [33]). Actually, all existing neutrino oscillation data, except the 
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17.6% 
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Table 2: The 3a allowed ranges and % spread of sin 2 # 12 obtained in a two-neutrino oscillation 
analysis of the solar neutrino data, including the simulated data on the pp neutrino flux, and the 
KamLAND 766.3 Ty spectrum data. 

data of LSND experiment 9 [34], can be described assuming 3- neutrino mixing. This warrants a 
3-neutrino oscillation analysis of the potential sensitivity of a LowNu pp neutrino experiment to 
sin 2 6*i2- 

The 3-neutrino oscillations of interest are characterized by the neutrino mass-squared dif- 
ferences which drive the solar and atmospheric neutrino oscillations, Am 2 x = Am^ > and 
Arajj = Am 2 tm respectively, and by the 3 mixing angles in the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) neutrino mixing matrix, 6> 12 , 6*23 and 13 . In the standard parametrization of the 
PMNS matrix (see, e.g., [33]), the angles 6>i 2 and 6> 23 coincide with the angles which control the 
oscillations of solar and atmospheric neutrinos, while #13 is the angle limited by the data from the 
CHOOZ and Palo Verde experiments [36]. The precise limit on #13 depends strongly on Am 2 im 
(see e.g. [37]) The existing atmospheric and reactor neutrino data imply [12, 13] 

sin 2 13 < 0.05, 99.73% C.L. (12) 

The aim of the analysis which follows is to quantify the uncertainty which the absence of precise 
knowledge of the value of the CHOOZ angle #13 introduces in the precision of sin 2 12 determination. 

The analyses of the latest SK atmospheric neutrino and of the global solar neutrino data show 
also that [6, 12, 13] (Am 2 ^ = |Am 2 tm | ~ 2.1 x 10~ 3 eV 2 and Am 2 21 = Am% ~ 8 x 10~ 5 eV 2 . Thus, 
we have Ara^ << jAmlJ. Under this condition the probabilities of survival of the solar v e and 
of the reactor P e , relevant for the 3-neutrino oscillation interpretation of the solar neutrino and 

9 In the LSND experiment indications for oscillations — > v e with (Ato 2 )lsnd — 1 eV 2 were obtained. The 
LSND results are being tested in the MiniBooNE experiment [35]. 
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Figure 5: Comparison of A% 2 vs. sin 2 9 12 obtained from a two generation analysis (solid black 
line) and three generation analysis (dashed red line) analyses of the world neutrino data including 
the pp rate. The three panels show the results for the three illustrative values of the pp rate. We 
present the results for 1% error in the pp rate. For the three- neutrino oscillation analysis sin 2 9 13 
is allowed to vary freely. 



KamLAND data, have the form: 

Pi: S cos 4 £ 13 P e 2 ; + sin 4 £ 13 (13) 

where P 2 ^ is the corresponding probability of survival of v e or v e in the case of 2-neutrino mixing 
(see, e.g., [38]). For the reactor v e detected at KamLAND , P 2 ^ is given by eq. (4). For solar 
neutrinos, P 2 ^ = P 2 e g is the v e survival probability in the case of 2-neutrino oscillation [39, 31] 
in which the solar electron number density N e is replaced by iV e cos 2 #13. From eqs. (2), (5), (12) 
and (13) we get: 

P?:( 8 B) = cos 4 # 13 sin 2 # 12 , (14) 
P?:(pp) = cos 4 Ml-^sin 2 2M. (15) 

Note that for a given value of P^(pp), a non-zero value of #13 decreases the measured value of 
#12, for the low energy pp flux. On the other hand, for a given value of P^( 8 B ), a non-zero #13 
increases the measured value of #12, for the higher energy 8 B flux. 
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Figure 6: The pp rate as a function of sin 2 # 12 . The black solid line is the rate for sin 2 # 13 = 
while the red dot-dashed line is sin 2 # 13 = 0.05. We also show our illustrative pp rates of 0.68, 0.72 
and 0.77 on the plot. 



For Amjj << | Am 2 ! |, the probability relevant for the interpretation of data from the CHOOZ 
experiment is given by 

P?: C hooz « 1 " sin 2 20 13 sin 2 (Am 2 1 L/4E). (16) 

Note that the probability Peecnooz depends on Am| 1; unlike the probabilities relevant for the 
interpretation of the solar and KamLAND data. In the analysis which we have performed jAm^j 
was allowed to vary freely within the 3a allowed range given in [6] 10 . 

In Fig. 5 we show the dependence of A% 2 on sin 2 #12, obtained in a 3- neutrino oscillation 
analysis of the combined KamLAND, CHOOZ and solar neutrino data, including the simulated 
data on the pp neutrino flux. The results shown are for the three illustrative central values of 
R pp considered earlier, and for the case of 1% error in R pp . Except for sin 2 ^, all the other 
parameters, including sin 2 # 13 , were allowed to vary freely in the analysis. For comparison, results 
for sin 2 #i 3 = are shown in the same figure. 

We find that 

• For a relatively low value of the measured pp rate, R pp = 0.68, a non-zero sin 2 6'i 3 leads to 
smaller minimal and maximal allowed values of sin 2 9 12 ■ 

• For R pp = 0.72, the minimal allowed value of sin 2 9 12 diminishes, while the maximal allowed 
value remains unaffected. 

10 Details of our 3-neutrino oscillation analysis can be found in [40]. 
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Figure 7: Same as Figure 3 but for a three generation analysis including all neutrino data and 
keeping sin 2 #13 free. We also show the 3a allowed band of sin 2 9 12 from the current data. 

• If R pp = 0.77, both values are practically unaffected. 

To help explain these features, we plot in Fig. 6 the pp rate (cf. Eqs. (6) and (15)) as a 
function of sin 2 9\2 for sin 2 #13 = and sin 2 #13 = 0.05. The three horizontal lines in the figure 
correspond to the three central values of R pp considered in our analysis. We see from the figure 
that, 

1. For a given value of sin 2 # 12 , a non-zero sin 2 # 13 always reduces the predicted pp rate, 

2. For a given measured pp rate, a non-zero sin 2 #13 would reduce the measured value of sin 2 9\2, 

3. For a given measured pp rate there is always a limiting value of sin 2 (9 12 , such that for sin 2 9 12 
exceeding this value a non-zero sin 2 # 13 does not affect the allowed values of sin 2 6 I 12 . 

Points (1) and (2) imply, in particular, that if for a given value of sin 2 9 12 the pp rate predicted 
assuming 2-neutrino oscillations is larger than the measured pp rate, a non-zero $13 can improve 
the quality of the fit. Thus, values of sin 2 9u smaller than the minimal allowed one determined in 
a 2-neutrino oscillation analysis, could become allowed in the case of 3-neutrino oscillations due 
to non-zero sin 2 9 13 . As a consequence of point 3 we can conclude that if for a particular value 
of sin 2 #12, the pp rate predicted in the case of 2-neutrino oscillations is larger than the measured 
pp rate, a sin 2 9 13 which differs from substantially and further lowers the pp rate, would not be 
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Figure 8: Same as Figure 7 but with the projected error from the third phase of the SNO experi- 
ment included. We show the 3<r allowed band from current data. 



favored by the data. We use these features to explain Fig. 5. However it is to be borne in mind 
that Fig. 6 contains only the pp rates while in Figure 5 we show the results from a combined 
analysis of global solar and KamLAND data. 

• We see from Fig. 6 that the pp rate, predicted in the case of 2-neutrino oscillations, exceeds 
0.72 for sin 2 # 12 ^ 0.27. Therefore in this case a non-zero sin 2 ^ improves the quality of 
the fit for all values of sin 2 #12 ^ 0.27. Thus, relatively small values of sin 2 #12 which are 
disfavored by the 2-neutrino oscillation analysis, could become allowed if sin 2 #13 7^ 0. This 
would lead to a smaller minimal allowed value of sin 2 #12. For sin 2 Qyi > 0.27, the predicted 
2-neutrino oscillation rate is already lower than 0.72, and therefore a non-zero #13 would not 
change the maximal allowed value of sin 2 # 12 • 

• The predicted R pp for 2-neutrino oscillations can be greater than R pp = 0.77 only for rela- 
tively small values of sin 2 # 12 , which are strongly disfavored (if not ruled out) by the current 
data. Consequently, a non-zero sin 2 6*13 is not expected to make any impact on the sin 2 # 12 
determination if the measured R pp = 0.77. 

• In the case of a measured R pp = 0.68, the predicted 2-neutrino oscillation pp rate is larger 
than 0.68 for sin 2 # 12 ^ 0.34 and a non-zero value of sin 2 # 13 can improve the quality of the 
fit for a large range of values of sin 2 dyi . This explains why the minimal allowed value of 
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sin 2 9 12 diminishes considerably for sin 2 #13 7^ in Fig. 5. The maximal allowed value of 
sin 2 9 12 is also seen to reduce. The reason for this can be understood if one notes that the 
values of sin 2 # i2 favoured by the pp experiment, for sin 2 9 13 = 0, would correspond to a 
relatively higher value of Rcc I 'Rnc compared to that measured at SNO. If #13 is non-zero, 
then the same R pp could be produced at a lower value of #12. Since the predicted Rcc /Rnc 
in SNO is given by Eq. (14), this lower sin 2 9 12 coupled with non-zero #13 would help reduce 
the 8 B probability and the pp "data" could be "reconciled" with the SNO CC/NC data. 
Therefore for R pp = 0.68, the best-fit from a three generation analysis comes at a non-zero 
value of sin 2 9 13 and a lower value of sin 2 9 12 . Note that the value of the global Xmin f° r the 
three generation case (#13 7^ 0) is lower than the two-generation case (#13 = 0). However as 
noted before, since the predicted R pp is already lower than 0.68 for sin 2 9 12 *S 0.34, the pp 
experiment would force #13 = for these high values of sin 2 #i 2 . Thus above sin 2 9 12 k, 0.34 
the x 2 remains the same irrespective of whether 9 13 was kept free (three-generation) or fixed 
at zero (two-generation). However, since the global Xmin was lower for the three-generation 
fit, the very high values of sin 2 #i 2 , which were allowed in the two-generation analysis get 
disfavored by the three-generation fit. 

In Fig. 7 we depict the sin 2 12 sensitivity of the world neutrino oscillation data, including the 
sample data on the pp rate, as a function of % error in the measured pp rate in the case when 
#13 is kept free. In Table 3 we give the corresponding 3<r ranges of allowed values and spread of 
sin 2 #12 (columns 3 and 4). A comparison of this figure with Fig. 3 and the Table 3 with Table 
2 demonstrates clearly all the specific features associated with the three values of the pp rate 
discussed above: 

For R pp = 0.77, #13 7^ hardly makes any difference. 

For Rpp = 0.72, the minimal allowed value of sin 2 # 12 diminishes as a consequence of #13 7^ 0, while 
the maximal allowed value is unaltered. 

For Rpp = 0.68, the inclusion of 13 7^ in the analysis leads to a reduction of both the minimal 
and maximal allowed values of sin 2 9 12 . 

For Rpp = 0.72 and 0.68 cases we also note that the minimal value of sin 2 #12 increases as the error 
in R P p increases and the SNO data begins to have a greater influence on the fit." Thus, the effect 
of the uncertainty due to sin 2 #13 on the precision of sin 2 # i2 measurement decreases as the error 
in the pp rate increases. 

Figure 8 shows the corresponding sensitivity plot of sin 2 # i2 from a 3-neutrino oscillation anal- 
ysis of the global neutrino oscillation data including both the (hypothetical) pp rate data and the 
prospective data (and errors) from phase-Ill of the SNO experiment. The corresponding allowed 
ranges and spread of sin 2 9 12 is given in columns 5 and 6 of Table 3. 

From the expression of the 3-neutrino oscillation probability, eq. (13), we see that the factor 
cos 4 #13 acts like a "normalization constant" . Since the current 3cr limit on this parameter is 
sin 2 #13 < 0.05, one would get a ~ 10% uncertainty in and would expect similar uncertainty 
to appear in the value of sin 2 #12 determined using the pp rate. The actual increase in the sin 2 #12 
uncertainty due to sin 2 #13 is smaller than ~ 10%, typically being ~ 3% for the plausible values 
of the pp rate (R pp = 0.72) we have considered. Even for the limiting case of R pp = 0.68, the 
maximal increase is by ~ 6% 
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The main reason for this is that when we include both "low" and "high" energy experiments 
in the global analysis there are two conflicting trends. While a non-zero value of sin 2 13 would 
have a tendency to lower the value of sin 2 # 12 determined from a 2-neutrino oscillation analysis 
of the data from the pp experiment, but it would also have a tendency to increase the value of 
sin 2 9 12 determined from the data of SNO and SK (i.e., S B neutrino) experiments. At the lower 
bound, in general, pp tries to shift the fit to non-zero sin 2 #13 and hence lower values of sin 2 6>i 2 but 
the global data including SNO prevents that. On the other hand, at the upper bound, in general, 
SNO can push the fit to sin 2 9 13 ^ and higher sin 2 9 12 but pp prefers to keep it at sin 2 13 = 0.0 
which corresponds to a lower maximal allowed value of sin 2 12 . 

The uncertainty in the value of sin 2 #13 leads to an error in sin 2 9 12 at the few percent level only 
when sin 2 9 12 is determined using data on the pp rate together with the global solar and reactor 
neutrino data. Nevertheless, the spread in the value of sin 2 12 remains well above 12.5% and 
typically exceeds 16% at 3a. As we will show in the next section, sin 2 9 12 could be measured with 
a considerably higher precision in a reactor neutrino experiment with a baseline tuned to SPMIN. 
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Table 3: The 3a allowed ranges and % spread of sin 2 # 12 , obtained from a 3-neutrino oscillation 
analysis of the global solar and reactor neutrino data, including the hypothetical data on the pp 
rate. 



4 Measuring #12 in a Reactor Experiment at SPMIN 

In this section we investigate the possibility of measuring the solar neutrino mixing parameter 
sin 2 9 12 in a reactor v e oscillation experiment, in which the baseline is chosen to correspond to 
a minimum of the v e survival probability (SPMIN) [16]. We will consider in what follows an 
experiment similar to KamLAND , but with a baseline tuned to the SPMIN. The condition of 
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SPMIN reads 



L mm « 1.24 £ ( M '^„ m (17) 
Amji /e^ 2 

For the "old" low-LMA best-fit value of Ara^ = 7.2 x 10~ 5 eV 2 , the baseline which allows the 
most accurate measurement of sin 2 6\2 was found to be 70 km [16]. For these Am^ and baseline 
the SPMIN appears at the prompt e+ energy of E vis = 3.2 MeV 11 , E vis 9* E-0.8 MeV. The latter 
corresponds to the maximum of the e + (event) spectrum in the absence of oscillations. Obviously, 
the energy E vis ~ 3.2 MeV, is the most relevant for the statistics of the experiment. We will 
show in this section that for the current global best-fit value of Am^ = 8.3 x 10~ 5 eV 2 , sin 2 9 12 
could be measured with an accuracy of ~ 2% (6%) at la (3<x) if the baseline chosen is L ~ 60 
km. For Am^ = 8.3 x 10~ 5 eV 2 and L ~ 60 km, the SPMIN appears at E vis = 3.2 MeV in the 
e + spectrum. We extend our earlier work [16, 41, 18] by investigating in detail the dependence of 
the precision of sin 2 # 12 measurement on the true value of Am^, the baseline, the statistics and 
on the systematic error of the experiment. We obtain results assuming 2-neutrino oscillations and 
compare them with the results of a 3-neutrino oscillation analysis. In the latter sin 2 #13 is allowed 
to vary freely within its currently allowed range. We discuss also the relevance of the geo-neutrino 
flux for the precision of sin 2 #12 measurement. 



4.1 Sensitivity to sin 2 #12 and the Baseline of the Experiment 

In Fig. 9 we show the sensitivity to sin 2 #i 2 , expected in a reactor experiment, as a function of 
the baseline L. We assume a total systematic uncertainty of 2% and consider statistics of 73 
GWkTy (given as a product of reactor power in GW and the exposure of the detector in kTy). 
The detector material composition is assumed to be the same as that of KamLAND and so the 
detector considered has the same number of target protons per kton as KamLAND. The total 
reactor power, the detector size and the exposure time are kept the same for all baselines. Thus, 
for longer baseline the number of events would decrease as ~ L~~ 2 . We assume that the true value 
of sin 2 6*12 = 0.27 and simulate the prospective observed positron spectrum in the detector for four 
different assumed true values of Am^. This figure is obtained for sin 2 #13 = 0. 
We define a y 2 function given by 

x 2 = J2(N? ata - N t i heory )(al)- 1 (Nf ata - Nf eory ) , (18) 

where N" (a = data, theory) is the number of events in the i th bin, <r 2 - is the covariant error 
matrix containing the statistical and systematic errors and the sum is over all bins. We use this 
X 2 to fit the simulated spectrum data and get the "measured" value of sin 2 # 12 , keeping Am^ free. 
We simulate the spectrum at each baseline and plot the range of values of sin 2 6 12 allowed by the 
simulated data as a function of the baseline. The baseline at which the band of allowed values 
of sin 2 9 12 is most narrow is the "ideal" baseline for the SPMIN reactor experiment. The figure 
confirms that this "ideal" baseline depends critically on the true value of Am^ (cf. eq. (17)). 

11 For the details of our reactor code and statistical analysis see refs. [16, 41, 18]. 
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Figure 9: Sensitivity plots showing the la, 1.64a, 2a, and 3a range of allowed values for sin 2 6* 12 as 
a function of the baseline L. The 4 panels are for 4 different true value of Am|. The true value 
of sin 2 9 12 is assumed to be 0.27 in all the cases. The Am^ is allowed to vary freely in the fit. 

The optimal baseline for the true value of Am^ = 8.3 (8.0) x 10~ 5 eV 2 is seen from Fig. 9 to 
be 60 (63) km, while for the "old" low-LMA best-fit value of Am^ = 7.2 x 10~ 5 eV 2 the best 
baseline would be 70 km. At the optimal baseline the SPMIN reactor experiment can achieve an 
unprecedented accuracy of ~ 2% (6%) at la (3a) in the measurement of sin 2 #12 • 

Figure 9 suggests that the optimal baseline for a given true value of Am 21 is very finely tuned. 
For instance, if for Am 21 (tr"ue) = 8.3 x 10 -5 eV 2 we change the baseline from L = 60 to L = 50 
km, the sensitivity in sin 2 # 12 decreases from ~ 2% (6%) to ~ 3% (11%) at la (3a). However, note 
that Fig. 9 was obtained by allowing Am 21 to vary freely. This is equivalent to assuming that 
both sin 2 # 12 and Am^ are determined in the reactor SPMIN experiment. For some baselines, 
especially at smaller L, the oscillation induced spectral distortion is not large enough to measure 
Am 21 sufficiently accurately, while for the longer baselines the statistics is lower. These factors 
lead to a certain uncertainty in the determination of Am 21 with the experimental set-up under 
discussion. The uncertainty in the Am^ determination translates into additional uncertainty 
in the measured sin 2 # 12 . If Am 21 could be measured with a sufficiently high precision in an 
independent experiment, the uncertainty in sin 2 9\% due to Am 21 would be reduced. 

Figure 10 represents a sensitivity plot similar to that shown in Fig. 9, but obtained for Am 21 
fixed at its assumed true value (indicated on each of the panels). As Fig. 10 shows, for fixed Am 21 
assumed to have been determined with a sufficiently high precision in an independent experiment, 
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Figure 10: The same as in Fig. 9, but for fixed Am^ having an assumed true value indicated in 
each of the panels of the figure. 

the choice of the baseline for setting up the SPMIN experiment becomes broader. It follows from 
Fig. 10 that for Am^trwe) = 8.3 x 10~ 5 eV 2 , for instance, the change of the baseline from L = 60 
to L = 70 km, leads to a minor increase of the uncertainty in the value of sin 2 9\2 from 6.1% to 
6.3% at 3d. 

It is actually quite possible that Am^ will be measured with a rather high accuracy in the 
future. The KamLAND experiment could determine Am^ with an error of about 7% (at 3cr) 
using data of 3 kTy [28, 16, 12]. The suggested SK-Gd experiment [42] has the potential of 
measuring the value of Am 2 x with an error of ~ 2 — 3% (at 3a) [18]. In Fig. 11 we show the 
sensitivity to sin 2 9i2 expected if we combine the SPMIN reactor data with 3 kTy prospective 
data from KamLAND (lower left panel) and simulated 5 year data from the suggested SK-Gd 
experiment (lower right panel). The upper panels were obtained using data from the SPMIN 
reactor experiment alone. For all the panels we have assumed Am 2 1 (trwe) = 8.3 x 10~ 5 eV 2 . In 
the upper left panel we allow Am^ to vary freely, while in the upper right panel Am^ is fixed at 
the assumed true value. We note that if the SPMIN reactor data is combined with 5 year data 
from the SK-Gd experiment, the choice of optimal baseline is much wider since Am^ would be 
determined with a relatively high precision by the SK-Gd experiment. With the addition of the 
SK-Gd results to the total data set, the spread in sin 2 # 12 is ~ 5.7% at 3a. The combined SPMIN 
reactor and KamLAND 3 kTy data would yield an uncertainty in the value of sin 2 # 12 of ~ 5.9% at 
3a. Since the analysis of the combined KamLAND (or SK-Gd) and SPMIN reactor data confirms 
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Figure 11: Sensitivity plots showing the la, 1.64a, la, and 3a range of allowed values for sin 2 # 12 as 
a function of the baseline L. The true values of Am^itrue) = 8.3 x 10 -5 eV 2 and of sin 2 #12 (true) = 
0.27 for all panels. The upper left hand (right hand) panel is obtained allowing Am^ to vary 
freely (for fixed Am^trwe) = 8.3 x 10~ 5 eV 2 ). The results shown in the lower left hand (lower 
right hand) panel are from combined analysis of 3 kTy data from KamLAND (5 year data from 
the suggested SK-Gd experiment) and 73 GWkTy data from a SPMIN reactor experiment. 

that the effect of Am| x on the sin 2 #12 sensitivity can be negligible, we will take Am^ to be fixed 
for the remainder of this section. 

4.2 Impact of Statistical and Systematic Errors on sin 2 9u Sensitivity 

One of the important requirements for the type of high precision experiment we are discussing is 
the accumulation of relatively high statistics in a reasonable period of time. Since the statistics 
falls as ~ L~ 2 and since rather long baselines are required for a precision measurement of the solar 
neutrino oscillation parameters, for a given reactor power longer baselines would imply bigger 
detectors and larger exposure times. Thus, the question about the dependence of the precision 
of measurement of sin 2 #12 in a reactor SPMIN experiment on the statistics of the experiment 
naturally arises. In Fig. 12 we show the effect of the statistics on the sin 2 612 sensitivity in the 
case of Am 2 1 (trwe) = 8.3 x 10~ 5 eV 2 . The four panels are obtained for four different sample 
baselines of 50, 60, 70 and 80 km. The range of allowed values of sin 2 # 12 is shown as a function 
of the product of reactor power and the detector mass and exposure time. For L = 60 km, 
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Figure 12: Sensitivity plots showing the la, 1.64a, 2a, and 3a ranges of allowed values for sin 2 6 l2 
as a function of the statistics in units of GWkTy. All the four panels correspond to a fixed value 
of Am 2 21 = 8.3 x 10~ 5 eV 2 . 

for instance, the uncertainty in sin 2 9\2 diminishes from 3% (10%) to 2% (6%) at la (3a) as the 
statistics is increased from 20 GWkTy to 60 GWkTy. Note that the difference in the sin 2 # 12 
precision for 60 GWkTy and 73 GWkTy (used in Figs. 9 and 10) is marginal, and shows up only 
in the first place in decimal in the value of the spread. 

Another important aspect which determines the potential of the experiment for precision mea- 
surement of sin 2 9 12 is the systematic uncertainty. Obviously, smaller systematic errors are prefer- 
able. All the plots presented so far in this section have been generated with an assumed 2% 
systematic error. The systematic uncertainty in the KamLAND experiment is about 6.5%. Most 
of it comes from the uncertainty in the detector fiducial mass and the reactor power. Our choice 
of 2% for the systematic error is based on the optimistic assumption that the error in the flux nor- 
malization could be reduced sufficiently by using the near-far detector set-up. One could envisage 
the 9 12 reactor SPMIN experiment as a second leg of a reactor experiment dedicated to measure 
9i 3 (see, e.g., [43]). The detector for measuring # 13 could then effectively be used as near detector 
for the long baseline SPMIN experiment for high precision measurement of 6 , 12 . It should be added 
that the errors due to the uncertainties in the threshold energy and v e spectrum have also to be 
reduced to achieve the systematic error of 2%. Experimentally this could be a very challenging 
task. 

Since systematic uncertainties may be difficult to reduce in the experiment under discussion, 
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Figure 13: The Same as in Fig. 10 but for a systematic uncertainty of 5%. 

we estimate next how much the precision on sin 2 9\2 deteriorates as the systematic error increases. 
Figure 13 shows the effect of increasing the systematic error from the 2% assumed by us to the 
rather conservative value of 5%. For Am^itrue) = 8.3 x 10~ 5 eV 2 , the spread in sin 2 # 12 at 
L = 60 km increases from 6.1% to 8.6% at 3a, as the systematic error is increased from 2% to 5%. 
We conclude that the effect of systematic uncertainty on the precision of sin 2 612 measurement 
is important, but its impact is not dramatic as long as the systematic error does not exceed 
5%. Similar conclusion regarding the effect of a 4% systematic error on the accuracy of sin 2 #12 
determination in the SADO experiment was reached in [23]. 

4.3 The Uncertainty Due to sin 2 #13 

As we have discussed earlier in connection with the KamLAND experiment, the 3- neutrino oscil- 
lation survival probability for the reactor v e of interest is given by 

P ee « cos 4 13 ^1 - sin 2 26 12 sin 2 ^ , (19) 

where the term ~ sin 4 #13 has been neglected. Therefore the uncertainty in sin 2 #13, eq. (12), 
brings up to a ~ 10% uncertainty in the value of the v e survival probability. Since the factor 
cos 4 #i3 can only reduce the survival probability, it does not affect the upper limit of the allowed 
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Figure 14: The same as in Fig. 10 but for a 3-neutrino oscillation analysis in which sin 2 13 is 
allowed to vary freely within its current 3a allowed range. 



range of sin 2 12 . However, it can have an effect on the minimal allowed value of sin 2 i2 reducing 
it further, and thus can worsen, in principle, the precision of the experiment to sin 2 12 . 

The additional error on sin 2 20i2 coming from the uncertainty in sin 2 #13 can be roughly esti- 
mated using eq. (19) as [41], 
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where AP ee and A(sin 2 0i3) are the uncertainties in the determination of the survival probability 
and sin 2 0i3, respectively. In the SPMIN region we are interested in one has sin 2 (Am 21 L/4i? i/ ) ~ 1 
and therefore 



5(sin 2 20i 2 ) « 2AP ee sin 2 13 + 2 cos 2 29 12 A(sin 2 6 13 ) . (21) 

Thus, for a reactor SPMIN set-up, the first term gives an extra contribution of about 2AP ee sin 2 #13 
to the allowed range of sin 2 20 12 . Even under the most conservative conditions one can expect that 
AP ee ^ 0.1, so this term could give an additional contribution of ^ 0.01 to the allowed range of 
sin 2 2012. The second term is independent of the precision of a given experiment. It depends only 
on the best-fit value of cos 2 20i2 and on the uncertainty in sin 2 0i3. For the current 3a upper limit 
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on sin 2 9 13 of 0.05 and best-fit value of cos 2 29\ 2 = 0.19, the second term would lead to an increase 
in the uncertainty in sin 2 29 12 by about 0.02 only. The suppression of this term is mainly due to 
the presence of the cos 2 29i 2 factor, which is relatively small for the current best-fit value. Thus, 
even though the current uncertainty in sin 2 # 13 brings a 10% uncertainty in the value of P ee , it 
increases the allowed range of sin 2 29± 2 only by a few %, if one uses a reactor experiment with a 
baseline tuned to the SPMIN for a high precision measurement of sin 2 29\ 2 . 

The above conclusions are illustrated in Fig. 14 showing the uncertainty in sin 2 9 12 expected in 
the case when 9 13 is allowed to vary freely in its currently allowed range of sin 2 9 13 < 0.05. The fig- 
ure confirms that the upper bound on sin 2 9\ 2 remains unaffected by the sin 2 # 13 uncertainty, while 
the minimal allowed value diminishes somewhat, increasing the uncertainty in sin 2 #12. However, 
for the baseline which corresponds to the SPMIN, the sensitivity reduces only by 2 — 3% in spite 
of the 10% uncertainty in sin 2 #13. For Aml^true) = 8.3 x 10 -5 eV 2 , for instance, the uncertainty 
in sin 2 9\ 2 increases from 6.1% to 8.7% at 3a. 

4.4 On the Impact of Geo-Neutrino Flux 

Our Earth is known to be a huge heat reservoir and is estimated to radiate about 40 TW of heat. 
A large fraction (~ 16%) of this is believed to be radiogenic in origin, coming from the decay 
chain of 238 U, 232 Th and 40 K. The radioactive decays of these isotopes produce antineutrinos in 
the beta decay processes of their decay chains. These v e coming from inside the Earth are usually 
called Geo-neutrinos {v 9eo ) [44]. The maximum energy of the v 9 e eo produced in the 40 K decay 
chain is only E Pe = 1.31 MeV, which is below the detection energy threshold of v e in scintillation 
detectors of the type of KamLAND we are considering. However, the u 9eo from 238 U and 232 Th 
have maximum energy of E Pe = 3.26 MeV and E Pe = 2.25 MeV, respectively, and can be observed 
in scintillation detectors. 

The flux of v 9 e eo is unknown. Even the total heat radiated by Earth has a rather large uncer- 
tainty: it could be 31 — 40 TW. There is no direct measurement of the abundances of the 238 U 
and 232 Th inside the Earth. One can estimate their abundances using the meteoritic and seismic 
data. This results in the v 9 e eo flux being largely model dependent and uncertain. Most models 
give the bulk 232 Th/ 238 U ratio as 232 Th/ 238 U ~ 3.8. Even the value of this ratio could have a 
realtively large error (e.g., the authors of [45] estimate this error as 14%). The measurement of the 
P 9eo flux would lead to a better understanding of the interior of the Earth, and is therefore a very 
important branch of neutrino physics in its own right [46]. As far as the precision measurement 
of the neutrino oscillation parameters is concerned, the events due to v 9 e eo can be an important 
background and can lead to an error in the measured value of sin 2 9\ 2 . 

Since the v 9 e eo have a maximum energy of E De = 3.26 MeV which corresponds to a prompt e + 
energy of only E vis = 2.48 MeV, one way to avoid the uncertainty due to v 9eo is to implement a 
prompt energy threshold of 2.6 MeV, as is done by the KamLAND collaboration. In this paper 
we have followed the KamLAND approach. In this case the observed e + spectrum does not have 
any "contamination" due to contributions from v 9eo . An alternative approach is to use the entire 
prompt e + energy spectrum in the analysis, taking the i> 9eo flux into account. Since the theoretical 
estimates on the u 9eo flux are presently rather imprecise, one could let the 238 U and 232 Th v 9eo 
flux normalisation vary as a free parameter. Both approaches have their merits and drawbacks. 
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In the first approach (we use in this paper), while there are no additional uncertainties due to the 
unknown u 9eo background, one has to contend with the experimental challenge of understanding 
and reducing the error associated with the prompt e + threshold energy. In the case of KamLAND 
experiment the uncertainty in the e + threshold energy corresponds to a systematic error ~2.3%. In 
the second approach there is no prompt e + energy cut, but one has to handle the uncertainty due to 
the lack of knowledge of v 9 e eo background. Keeping the 238 U— and 232 Th— v 9 e eo flux normalisation 
as free parameter brings in extra error in the measurement on sin 2 6> 12 . 

The key feature in the 9 12 SPMIN reactor experiment suggested in [16], is the appearance of 
SPMIN in the observed e + spectrum. If the SPMIN appears at a prompt energy of E vis > 2.6 
MeV, implementing a threshold of E V i S = 2.6 MeV 12 and thus avoiding the v 9 e eo background 
might permit to measure sin 2 d\2 with the highest precision, achievable in the experiment under 
discussion. If, on the other hand, SPMIN appears at E vis < 2.6 MeV, the entire e + energy 
spectrum would have to be taken into account and in this case the v 9 e eo background cannot be 
avoided. For a given value of Am^, the position of the SPMIN in the v e spectrum depends on the 
baseline of the experiment. For shorter baselines, SPMIN occurs at smaller energies. Therefore 
the choice of the baseline of the experiment would determine whether one would have to take the 
v 9eo background into account or not. 

The authors of [23] have included the u 9eo background in their analysis of the sin 2 Qyi precision 
expected in the SADO experiment in Japan with a baseline of L = 54 km. They conclude that for 
this experimental set-up, the v 9 e eo background does not have significant impact on the precision 
of sin 2 #12 measurement. For L = 54 km, the SPMIN is at E vis = 2.8 MeV. The uncertainty in 
the z/| eo flux does not make much impact in this case since E vis = 2.8 MeV is larger than the 
background v 9eo energies. For shorter baselines the SPMIN will take place at E vis < 2.6 MeV and 
the uncertainty due to the u 9eo flux can affect noticeably the precision of measurement of sin 2 Qyi- 

5 Conclusions 

We have investigated the possibilities of high precision measurement of the solar neutrino mix- 
ing angle #12 in solar and reactor neutrino experiments. As a first step, we have analyzed the 
improvements in the determination of sin 2 #i 2 , which can be achieved with the expected increase 
of statistics and reduction of systematic errors in the currently operating solar and KamLAND 
experiments. With the phase-Ill prospective data from SNO experiment included in the current 
global solar neutrino and KamLAND data, the uncertainty in the value of sin 2 #12 is expected to 
diminish from 24% to 21% at 3a. If instead of 766.3 Ty, one uses simulated 3 kTy KamLAND 
data in the same analysis, the 3a error in sin 2 6*12 reduces to 18%. 

We next considered the potential of a generic LowNu v — e elastic scattering experiment, 
designed to measure the pp solar neutrino flux, for high precision determination of sin 2 # 12 . We 
examined the effect of including values of the pp neutrino induced electron scattering rates in 
the x 2 analysis of the global solar neutrino data. Three representative values of the rates from 

12 This will increase the systematic uncertainty due to the error from the prompt e + energy cut. We have shown 
that the impact of the increase of the systematic uncertainty on the precision of sin 2 6*12 determination is relatively 
small as long as the systematic error does not exceed 5%. 
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the currently allowed 3a range were considered: 0.68, 0.72, 0.77. The error in the measured 
rate was varied from 1% to 5%. By adding the pp flux data in the analysis, the error in sin 2 9 12 
determination reduces to 14% (19%) at 3a for 1% (3%) uncertainty in the measured pp rate. 
Performing a similar three-neutrino oscillation analysis we found that, as a consequence of the 
uncertainty on sin 2 #13, the error on the value of sin 2 9i 2 increases correspondingly to 17% (21%). 

We also studied the possibility of a high precision determination of sin 2 #12 in a reactor exper- 
iment with a baseline corresponding to a Survival Probability MINimum (SPMIN). We showed 
that in a L ~ 60 km experiment with statistics of ~60 GWkTy and systematic error of 2%, sin 2 9 12 
could be measured with an uncertainty of 2% (6%) at la (3a). The inclusion of the sin 2 6* 13 uncer- 
tainty in the analysis changes this error to 3% (9%). An independent determination of Am^ with 
sufficiently high accuracy would allow, as we have shown, sin 2 612 to be measured with the highest 
precision over a relatively wide range of baselines. We investigated in detail the dependence of 
the precision on sin 2 612 which can be achieved in such an experiment on the baseline, statistics 
and systematic error. More specifically, with the increase of the statistics from 20 GWkTy to 60 
GWkTy, the error diminishes from 3% (10%) to 2% (6%) at la (3a). For statistics of (60 - 70) 
GWkTy, the increase of the systematic error from 2% to 5% leads to an increase in the uncertainty 
in sin 2 12 from 6% to 9% at 3a. 

We have found that the effect of sin 2 # 13 uncertainty on the sin 2 6>i 2 determination in LowNu 
pp and SPMIN reactor experiments considered is considerably smaller than naively expected. 

The results of our analyses for the currently running, the proposed LowNu and future reactor 
experiments show that the most precise determination of sin 2 d\ 2 can be achieved in a dedicated 
reactor experiment with a baseline tuned to SPMIN associated with Am^ = Am|. 
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